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Abstract 

This contribution discusses the phenomenon of parton saturation, the 
color glass picture of hadronic wavefuntions, and their relevance in the 
early stages of nucleus- nucleus collisions. Evidence for the color glass 
condensate in the presently available RHIC data is critically reviewed. 
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1 Introduction 



The degrees of freedom involved in the early stages of a nucleus-nucleus collision 
at sufficiently high energy are partons, mostly gluons, whose density grows as 
the energy increases (i.e., when x, their momentum fraction, decreases). This 
growth of the number of gluons in the hadronic wave functions is a phenomenon 
which has been well established at HERA. One expects however that it should 
eventually "saturate" when non linear QCD effects start to play a role. 

The existence of such a saturation regime has been predicted long ago, to- 
gether with estimates for the typical transverse momenta where it sets in. But 
it is only during the last decade that equations providing a complete dynamical 
description of the saturated regime have been obtained. A remarkable feature 
which emerges from the solution of these equations is that the dense, saturated 
system of partons to be found in hadronic wave functions at high energy has 
universal properties, the same for all hadrons or nuclei. It follows that the 
early stages of hadronic collisions at sufficiently high energies are governed by 
universal wave functions ( "color glass condensate" ) , whose properties could in 
principle be calculated from QCD. This is a very exciting perspective which 
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fully justifies the active search for evidence of this novel regime of QCD both 
at HERA and at RHIC. 

It is expected that the saturation regime sets in earlier (i.e., at lower energy) 
in collisions involving large nuclei than in those involving protons. In fact, the 
parton densities in the central rapidity region of a Au-Au collision at RHIC 
are not too different from those measured in deep inelastic scattering at HERA. 
There is however one important difference: while at HERA these densities result 
from gluon evolution, at RHIC there is little evolution, at least in the central 
rapidity region, and the large densities result mostly from the additive contri- 
butions of the participant nucleons. Of course, one has the possibility at RHIC 
to explore various situations. In particular the study of dA collisions in the 
fragmentation region of the deuteron gives access to a regime where final state 
interactions should play a minor role and where quantum evolution could be 
significant. Indeed, very exciting results have been obtained in this regime. 

One may worry that at present energies, the conditions for saturation are 
at best marginally satisfied. Nevertheless a successful phenomenology based 
on the saturation picture has been developed over the last few years, both at 
RHIC and at HERA. We shall review in this report some of the experimental 
findings providing evidence for the color glass picture. It should be emphasized 
that the field is still in a rapidly evolving phase, most theoretical analysis are 
incomplete, so that conclusions reached today can at best be considered as 
tentative. Before we go into this phenomenological discussion, we shall start 
by a brief (and incomplete) historical perspective, to emphasize the remarkable 
theoretical developments of the last decade. 

2 The infancy of the idea of saturation 

An important feature of partonic interactions is that they involve mostly par- 
tons with comparable rapidities. Let us then consider a nucleus- nucleus collision 
in the target rest frame and consider what happens when one boosts the pro- 
jectile, increasing its rapidity in successive steps. In the first step, the valence 
constituents become Lorentz contracted in the longitudinal direction while the 
time scale of their internal motions is Lorentz dilated. In addition, the boost 
reveals new vacuum fluctuations coupled to the boosted valence partons. Such 
fluctuations are not Lorentz contracted in the longitudinal direction, and repre- 
sents the dynamical degrees of freedom. Making an additional step in rapidity 
would freeze these fluctuations, while making them Lorentz contracted as well. 
But the additional boost also produces new quantum fluctuations, which be- 
come the new dynamical variables. This argument can be repeated, and one 
arrives at the picture of a high-energy projectile containing a large number of 
frozen, Lorentz contracted partons (the valence partons, plus all the quantum 
fluctuations produced in the previous boosts), and partons which have a small 
rapidity and are not Lorentz contracted. This space-time description was de- 
veloped before the advent of QCD (see for instance [1]; in Bjorken's lectures 
[2], one can actually foresee the modern interpretation of parton evolution as a 
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rcnormalization group evolution a la Wilson on which the color glass formalism 
is based). 

Of course, such a space-time picture, which was deduced from rather general 
field theoretical considerations, can now be understood in terms of QCD. In 
fact, shortly after QCD was established as the theory of strong interaction, 
quantitative equations were established, describing the phenomenon outlined 
above [3-8]. In particular, the equation derived by Balitsky, Fadin, Kuraev and 
Lipatov [3,4] describes the growth of the non- integrated gluon distribution in a 
hadron as it is boosted towards higher rapidities. Experimentally, an important 



Figure 1: The gluon structure function in a proton measured at HERA. 

increase of the number of gluons at small x has indeed been observed in the 
Deep Inelastic Scattering experiments performed at HERA (see fig. 1), down to 
x ~ 10~ 4 . Such a growth raises a problem: if it were to continue to arbitrarily 
small x, it would induce an increase of hadronic cross-sections as a power of the 
center of mass energy, in violation of unitarity bounds. 

However, as noticed by Gribov, Levin and Ryskin in [9], the BFKL equation 
includes only branching processes that increase the number of gluons (g — * gg 
for instance) , but not the recombination processes that could reduce the number 
of gluons (like gg — ► g). While it may be legitimate to neglect the recombination 
process when the gluon density is small, this cannot remain so at arbitrarily high 
density: a saturation mechanism of some kind must set in. Treating the partons 
as ordinary particles, one can get a crude estimate of the onset of saturation 
from a simple mean free path argument. The recombination cross-section for 
gluon with transverse momentum Q is roughly given by 



Hl+ZEUS 



» 20 




Q 2 =20 GeV 2 



HI NLO-QCD Fit 2000 



X 



a 



Q 2 



(1) 



3 



while the number of such gluons per unit of transverse area is given by 



P 



xG(x,Q 2 ) 



(2) 



where R is the radius of the hadron and x the momentum fraction of the con- 
sidered gluons. Saturation sets in when per <~ 1, or equivalently for: 



The momentum scale that characterizes this new regime, Q Sl is called the sat- 
uration momentum [10]. Partons with transverse momentum Q > Q s are in a 
dilute regime; those with Q < Q s arc in the saturated regime. 

A more refined argument for the onset of saturation was given in [11], where 
recombination is associated with a higher twist correction to the DGLAP equa- 
tion. More generally, one may view the saturation momentum Q s as the scale 
at which QCD non linear effects become important. This occurs when field 
fluctuations are such that gluon kinetic energies become comparable to their 
interaction energies, that is when ((dA) 2 ) <~ a s ((A 2 ) 2 ), where (A 2 ) denotes the 
fluctuations of the gauge fields. Since the relevant dynamics is in the transverse 
plane, the magnitude of the gradient is fixed by the transverse momentum Q. 
As for the magnitude of the gauge field fluctuations, (A 2 ), they can be esti- 
mated from the particle number density in the transverse plane, i.e. (A 2 ) <~ p, 
with p given in eq. (2). The condition ((dA) 2 ) <~ a s ((A 2 ) 2 ) translates then into 
Q 2 <~ a s (A 2 ) 1 which is eq. (3) above. Note that at saturation, naive perturba- 
tion theory breaks down, even though a s (Q s ) may be small if Q s is large: the 
saturation regime is a regime of weak coupling, but large density. 

Early estimates of Q s in nucleus-nucleus collisions were given in [12], and do 
not differ much from more modern ones [13]. But, as we have already empha- 
sized, what has changed dramatically over the last decade is that we do not have 
only access to the boundary of the saturation region, we now have a complete 
dynamical picture of the saturation region. Before we turn to a short introduc- 
tion to these new developments, let us indicate some expected characteristics of 
the saturated regime. 

Note first that the saturation momentum increases as the gluon density 
increases. This may come from an increase of the gluon structure function 
as x decreases. The increase of the density may also come from the coherent 
contributions of several nucleons in a nucleus. In large nuclei, one expects 
Q 2 oc ctgA 1 / 3 , where A is the number of nucleons in the nucleus. 

There is another feature. Consider the number of partons occupying a small 
disk of radius l/Q s in the transverse plane. This is easily estimated by com- 
bining eqs. (2) and (3); one finds this number to be proportional to l/a s . In 
such conditions of large numbers of quanta, classical field approximations may 
become relevant to describe the nuclear wave- functions. 



Q 2 = Q 2 S , with Q 2 S ~ a s (Q 2 s ) 



xG(x,Q 2 s ) 
ttR 2 



(3) 



4 



3 Modern formulation: 

the color glass condensate 

Once one enters the saturated regime the evolution of the parton distributions 
can no longer be described by a linear equation such as the BFKL equation. One 
of the major breakthrough of the last ten years is that non linear equations have 
been obtained which allow us to follow the evolution of the partonic systems 
form the dilute regime to the dense, saturated, regime. These take different, 
equivalent, forms, generically referred to as the JIMWLK equation. 

The color glass formalism, which provides the most complete physical pic- 
ture, relies on the separation of the degrees of freedom in a high-energy hadron 
into frozen partons and dynamical fields, as discussed above. In the original 
McLerran and Venugopalan model [14-16], the fast partons are frozen, Lorentz 
contracted, color sources flying along the light-cone, and constitute a density 
of color charge p(x±). Conversely, the low x partons are described by classi- 
cal gauge fields A M (x) determined by solving the Yang-Mills equations with the 
source given by the frozen partonic configuration. An average over all acceptable 
configurations must be performed. 

The weight of a given configuration is a functional W Xo [p] of the density p 
which depends on the separation scale xq between the modes which are described 
as frozen sources, and the modes which are described as dynamical fields. As 
one lowers this separation scale, more and more modes are included among the 
frozen sources, and therefore the functional W Xo evolves with xq according to a 
renormalization group equation [17-26]. 

This evolution equation for W^Jp] has been derived in [17-26] and reads 



where the kernel Xab(x±,y±) depends on the color density p only via Wilson 



Here V denotes an ordering along the x~ axis, A + is the classical color field of 
the hadron moving close to the speed of light in the +z direction. The field A + 
depends implicitly on the frozen sources, i.e. on the color charge density p(x±). 

This functional evolution equation can be rewritten as an infinite hierarchy 
of equations for correlation functions of the p's, or equivalently of the U's. For 
instance, the correlator tr(f/t (x±)U(y ± )) of two Wilson lines has an evolution 
equation that involves a correlator of four Wilson lines. If one assumes that this 
4-point correlator can be factored into a product of two 2-point functions, one 
obtains a closed equation for the 2-point function, called the Balitsky-Kovchegov 
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[19,21] equation: 



dtv(W(x ± )U(y ± )) a 



Us f 



(x± -y ± ) 2 



dln(l/x ) 2tt 2 J z± ( X± _ z ± ) 2 (y ± - z^f 

N c tr(rf(x ± )U(y ± )) xn tv(U\x ± )U(z ± ))tv(U\z ± )U( yi _)) 3 



(6) 



When the density of color charge p is small, one can expand the Wilson line 
U in powers of p. Eq. (6) becomes then 2 a linear evolution equation for the 
correlator (p(x±)p(y ± )} x , or equivalently for the unintegrated gluon density, 
which is nothing but the BFKL equation. 

The opposite limit, at very small Xq (where the color charge density is large - 
far inside the saturation region), also leads to significant simplifications. Indeed, 
since the exponent in the Wilson lines is then a large quantity, one can use a 
Random Phase Approximation in which the correlators of Wilson lines are small. 
The solution is a Gaussian (albeit non-local in the transverse coordinates) [27]: 



W X0 [p] = exp 



/ 



p{x ± )p(y A 



2p?(x±_ - y ± ,x ) 



(7) 



0, that is 



This formula is not valid to arbitrarily short distances \x± — y ± 
the domain of high Q 2 physics controlled by DGLAP evolution. 

Like with the BFKL or DGLAP evolution equations, the initial condition 
for the evolution is truly a non-perturbative input. One can in principle try to 
model it, and then adjust the parameters of the model to fit experimental data. 
A simple model is that proposed by McLerran and Venugopalan, in which the 
initial W xo [p] is a local Gaussian: 



W xo [p] = exp 



p(x±)p(x±) 



2p? 



(8) 



At this point, one should stress that testing the predictions of the Color Glass 
Condensate in principle requires to test both the properties of the evolution with 
rapidity and the initial condition. 

It is also important to note that the gaussians of eqs. (7) and (8) do not have 
the same status in the CGC framework: cq. (8) is one particular model for the 
initial condition at moderately small x, while eq. (7) is the asymptotic regime 
reached at very small x regardless of the initial condition. The latter is therefore 
a property of the small x evolution itself. The MV model requires an infrared 
cutoff at the scale A 



QCD ■ 



This is because assuming a truly local Gaussian 
distribution ignores the fact that color neutralization occurs on distance scales 
smaller than the nucleon size (~ A~* D ): two p's can only be uncorrelated if they 
are at transverse coordinates separated by at least the distance scale of color 
neutralization. In the asymptotic solution of cq. (7) , there is no infrared problem 
because the color neutralization is built in the dependence of p 2 (x± — y^,xo) 

2 The same is true of eq. (4), because in this limit the kernel \ab becomes quadratic in p. 
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on the transverse coordinates. Color neutralization in fact occurs on distance 
scales of the order of Qj 1 (xo) [27]. This is the physical origin of the universality 
of the saturated regime. 

4 Phenomenological predictions 

4.1 Geometrical scaling 

In the (ln(l/a;), ln(pj_)) plane, the saturation domain is defined by the condition 
^-qcd ^ P-L ~ Qs{x). In this region, one expects that observables (made 
dimensionless by scaling out the appropriate power of p±) depending a priori 
on both x and p± scale like functions of the unique variable p±/Q s (x). This 
property is called "geometrical scaling" in the literature. It has been shown 
by Iancu, Itakura and McLerran in [28] that this scaling property holds in a 
much larger domain, whose upper bound is p± < Q s (x) 2 /A QCD . This region 
has been labeled "extended scaling region" in figure 2. Technically, the reason 

Color Glass Extended 




ln(Pr) 

Figure 2: The scaling region in the (ln(l/a;),ln(px)) plane. 

why this scaling survives outside of the saturation domain is that the (linear) 
BFKL equation that controls the evolution in x when p± > Q s (x) preserves for 
a while the scaling properties inherited from the physics of the saturation region 
(which affects the BFKL evolution via a boundary condition at p± = Q s (x)). 

Such a scaling has been searched for in the data of the DIS experiments at 
HERA, and it turns out that one can indeed represent all the small x (x < 
10~ 2 ) data points for the structure function F 2 (:r,Q 2 ) as a single function of 
the variable r = Q 2 /Q 2 (x) with 

Q 2 s {x)^Q 2 (^) X . (9) 

This is illustrated by the plot of figure 3. The parameter Qo is set to 1 GeV, 
while xo w 3.10~ 4 and A « 0.29 are determined through a fit to the data [29]. 
By studying BFKL evolution to next-to-leading order, Triantafyllopoulos has 
obtained a value of A in good agreement with that observed [30] . 
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Figure 3: The -f*p cross-section at HERA, plotted against r = Q 2 /Q 2 (x). 



4.2 Structure functions and diffraction in DIS 

In an appropriate frame, and at leading order in a s , the cross section for the 
interaction of a virtual photon with a proton takes the following factorized form: 



a^ p (x,Q 2 ) = J dz J d 2 r ± \t[;(Q 2 7 z,r ± )\ 2 <T dipo i c (x,r ± ) . 



(10) 



In this formula, ip(Q 2 , z, r±) is the Fock component of the virtual photon light- 
cone wave function that corresponds to a qq dipole of size r±; it depends on 
the invariant mass Q 2 of the photon, on the transverse size r± of the qq dipole, 
and on the fraction z of the photon longitudinal momentum taken by the quark. 
The other factor in this formula, <7di P oic(^, r ±), is the total dipole-proton cross- 
section. It can be expressed in terms of a correlator of Wilson lines: 



<Jdi P o\c{x,r±) = J- Jd 2 X ± tr 



i-(u(x ± + ^)u\x ± - 2 



T>>] • <"> 



where the average is taken over the color field of the proton. 

Several models for this dipole cross-section have been used in the literature 
in order to fit HERA data. Golec-Biernat and Wusthoff have used a very simple 
parameterization [31,32] 



0-dipoie(a;, r±) = cr I - e 



(12) 



which has led to good results for describing the data at x < 10 2 and moderate 
Q 2 . In this formula, the scale Q s (x) was taken to be of the form given in 
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cq. (9). This model fails however at large Q 2 . This aspect was improved in 
[33], where the dipolc cross-section is parameterized in a way that reproduces 
pQCD for small dipoles. Note that these approaches, even if they are inspired by 
saturation physics, do not derive the dipole cross-section from first principles. 
Recently, Iancu, Itakura and Munier [34] derived an expression of the dipole 
cross-section from the color glass condensate framework, and obtained a good 
fit of HERA data with a small number of free parameters. An equally good fit 
has been obtained by Gotsman, Levin, Lublinsky and Maor who derived the x 
dependence of the dipole cross-section by solving numerically the BK equation 
and including DGLAP corrections [35]. 

One can express in a similarly factorizcd form the diffractive "/*p cross- 
section, where the exchange between the dipole and the target is color singlet. 
In fact, since the dipole cross-section involved in eqs. (10) and (11) is the total 
cross-section, one can use the optical theorem to deduce from it the forward 
elastic scattering amplitude: 

M elastic (x,r ± ,t = 0) = ^- J d 2 X ± tv [i-(u(X ± + ?±)U*(X ± -?±)) 

(13) 

assuming that this amplitude is predominantly imaginary. This can then be used 
to obtain for the diffractive DIS cross-section at t = an expression involving 
the square of the total dipole-proton cross-section 3 : 



da , Jx,Q 2 



dt 



^ d-'r ± \i,(Q\z,r ± )\ 2 ai ipole (x,Q*). (14) 

t=o 



By using in this formula the total dipole cross-section fitted to the inclusive 
F 2 , one obtains a good description of the measured diffractive DIS cross-section 
[31,32]. Note however that keeping only the qq component of the photon wave 
function is a good approximation only when the invariant mass of the diffrac- 
tively produced object is not too large. For large masses, one needs to consider 
also the qqg component. 



4.3 Particle multiplicities at RHIC 

In the framework of saturation models, predictions have been made for global 
observables like the total number of produced particles per unit of rapidity 
dN/drf. This can be studied as a function of the collision energy, the centrality 
of the collision, and the rapidity. 

The Color Glass Condensate may provide a dynamical justification rooted in 
QCD for many general features of particles production in hadronic interactions. 
This may be the case in particular for the phenomenon of "limiting fragmenta- 
tion", i.e., the expectation that the rapidity distribution in the fragmentation 

3 Another way in which this relation is stated in the literature is that in order to obtain 
the diffractive part of a cross-section, one should average the scattering amplitude over the 
configurations of the frozen sources before squaring the amplitude. 
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Figure 4: Particle multiplicity at RHIC for different collision energies, plotted 
against the rapidity of the produced particle relative to the rapidity of the beam 

rf = V - ?7bcam [36]. 

region becomes independent of the collision energy at high energy. Evidence for 
such a behavior has indeed been observed by the PHOBOS collaboration, as is 
illustrated in figure 4. Here, one has plotted the rapidity distribution dN/drf 
of produced particles in Au-Au collisions as a function of the rapidity relative 
to the incoming beam rapidity rf = r\ — ijbcam- When one increases the beam 
energy, one sees that dN/drf is the same at large rf for all energies. The par- 
ticles produced at some given (large) energy are those that would have been 
produced at lower energies plus new particles produced at a lower rapidity. An 
early interpretation of such a phenomenon in the language of the CGC has been 
proposed by Jalilian-Marian [37]. The interpretation suggested there is that in 
the fragmentation region of one of the nuclei (rf around zero) , this nucleus is a 
dilute partonic system while the other nucleus is in a saturated state. The rise 
in multiplicity when one decreases rf is then attributed to the growth of the 
parton distribution in the dilute nucleus. One should emphasize however that 
no detailed quantitative analysis of the data presented in figure 4 has yet been 
done in this framework. 

In comparing more quantitative predictions of the color glass formalism for 
dN/dr) to experimental data, it is important to keep the following points in 
mind. First, the color glass condensate can only predict the distribution of ini- 
tial gluons, set free typically at a proper time r ~ Qj - Between this early 
stage and the final chemical freeze-out, the system undergoes several non-trivial 
steps: kinetic and chemical equilibration (possibly with additional parton pro- 
duction), hadronization, etc., which are most often ignored in calculations based 
on the color glass condensate. Secondly, even the calculation of the initial gluon 
production is a highly non-trivial task. It involves, in principle, solving the 
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classical Yang-Mills equations in the presence of the color densities describing 
the distribution of frozen sources in both projectiles. This has been achieved 
analytically only for proton- nucleus collisions, when the color charge density 
inside the proton is assumed to be weak [38-40] . In the case of nucleus-nucleus 
collisions, two kinds of calculations have been performed: 

(i) Ab initio numerical calculations [41-44] that solve exactly the Yang-Mills 
equations for A A collisions. The average over the color charge densities in the 
projectiles is performed by assuming that the distributions W x „ [p] are given by 
the MV model, i.e. by the gaussian of eq. (8). Quantum evolution effects are 
therefore not included. Note that when one calculates dN/dn at rj = 0, the 
typical x probed in the wave function of the nuclei is about x ~ 1CP 2 . Such 
calculations can therefore be used to test the validity of the gaussian model 4 for 
W Xo [p] at xq ~ 10~ 2 as well as the non-linear dynamics of the classical field. 

(ii) Approximate analytical calculations of the initial gluon spectrum [45]. 
These calculations assume /cj_-factorization, although such a property is so far 
unprovcn for the collision of two dense projectiles. The non-integrated gluon 
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Figure 5: The number of particles per unit of rapidity in A A collisions, as 
a function of rapidity for different centralities (left plot) and as a function of 
centrality at different energies (right plot). From [45] . 

distribution in the nuclei is taken to be of the form <p(k±) ~ l/a s for k± <C 
Q s (x) and <p(k±) ~ Q^( x )/k± for k± » Q s (x) (saturation appears through the 
fact that the gluon distribution does not diverge like fcj 2 at small k±). The 
rapidity dependence in this model is governed by that of Q s (x), eq. (9), where 
the exponent A can be taken from the study of the scaling properties in the 
HERA data. The overall normalization factor, as well as the value of Q s at a 
certain fixed energy, were fitted on RHIC data at y/s = 130 GeV. A prediction 
was then made in [45] for the value of dN/drj at y/s — 200 GeV, which is in good 
agreement with experimental results, as illustrated in figure 5. One should also 
mention the fact that the value of Q 2 s (x) one needs in this approach in order 

4 Since the gluon multiplicity depends only on the correlator (A^ (x) A v (y)) , only one mo- 
ment of the distribution W XQ [p] is tested by this observable. One could certainly imagine 
different models for this distribution that give the same value for this particular correlator. 
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to reproduce the measured dN/dr], of the order of 2 GeV 2 , is relatively large 
compared to A QOD ; this may perhaps be taken as an encouraging indication for 
the validity of the overall picture. 

Note that the residual dependence on the number of participants in the 
quantity N~^. t dN/dr] calculated by Kharzeev, Levin and Nardi, comes entirely 
from the scale dependence of the strong coupling constant, as follows: 



1 dN 



Npart drj a s (Q 2 s ) 



(15) 



However, strictly speaking, since the gluon multiplicity is obtained by solving 
the classical Yang-Mills equations, there is no running of a s at this level of 
approximation. Certainly the running of a s will come together with next-to- 
lcading-order corrections, but it is at this point an ad hoc prescription. Since 
other interpretations of these same data, based on soft physics, are possible [46], 
it is unclear whether what we are "seeing" in the right hand panel of figure 5 
is indeed the running of a s {Q 2 s ) induced by the variation of the saturation scale 
with centrality. 
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Figure 6: The number of particles per unit of rapidity in dA collisions, as 
predicted in [47]. 

Such analysis of multiplicity distributions were extended to the case of 
deuteron-nucleus collisions [47] , with results in fair agreement with RHIC data 5 
as illustrated in figure 6. 



4.4 Spectra of produced particles 
4.4.1 p± dependence 

The saturation regime is characterized by a single scale, the saturation momen- 
tum Q s , and as was the case for DIS, one may look for scaling laws in various 

5 In the first version of this study, there was a discrepancy with the measured multiplicities, 
which seems to have been pinned down to a problem with the Glauber determination of the 
number of participants. 
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observables. An attempt has been made to identify such a scaling in the trans- 
verse mass distributions of produced particles [48]: it has been found that the 
spectra for different species of particles are well represented by a single function 
of m±/Q s . The centrality dependence of the scaling parameter Q 2 was found 
in rough agreement with Q 2 ~ ~ A 1 / 3 . It should be emphasized however, 

that the hadrons whose spectra are measured have certainly undergone many 
reinteractions, so that their spectra are unlikely to reflect directly the initial 
momentum distribution of the color glass. Thus the interpretation of the m± 
scaling in the saturation picture is unconvincing. 

In the case of deuteron-gold collisions, one does not expect final state in- 
teractions to play a dominant role. All the medium effects responsible for 
the difference between pA and pp collisions may in fact be taken into account 
in the color glass condensate. The produced partons hadronize then by the 
same mechanisms as in the vacuum, and their spectrum can be calculated by 
convoluting the partonic cross-section with the usual fragmentation functions. 
Moreover, since the calculation of the gluon spectrum only involves the corre- 
lator (U(0)U^(x±)) [39,49,50,40], one can in principle determine it from the 
dipole cross-section used in DIS [51]. A numerical calculation of the spectrum 
of hadrons produced in dA collisions has been performed in this approach in 
[52] , with results in fairly good agreement with the spectrum measured at RHIC 
in dA collisions. 

The situation is far more complicated in the case of nucleus- nucleus colli- 
sions. Indeed, the gluons which emerge from the color glass shortly after the 
initial impact (typically at a proper time of the order of Q^ 1 ~ 0.2 fm) will in- 
teract strongly and may form a hot and dense medium (a quark-gluon plasma?). 
These interactions will presumably modify the gluon spectrum: the momentum 
distribution will become more isotropic and the hard tail of the spectrum will 
be suppressed by parton energy loss. These additional effects are not taken into 
account in the color glass, which merely provides the initial condition for the 
subsequent evolution of the system. 

By assuming that local thermalization is achieved quickly, one can use hy- 
drodynamics to describe this evolution. One may adjust the initial conditions 
so that averaged quantities like the local energy density match those predicted 
in saturation models. Such a strategy has been used by Hirano and Nara [53] as 
well as Eskola, Niemi, Ruuskanen and Rasanen [54,55], who predicted hadron 
spectra in good agreement with the observed ones for p± < 2 GeV. One should 
however keep in mind that this kind of test is by construction only sensitive to 
integrated quantities predicted by the CGC rather than to the detailed shape 
of the spectrum: indeed, the very assumption that thermalization occurs means 
that the system has lost memory of the initial momentum distribution except 
for the local energy density. 

This remark implies that dA collisions are certainly much better suited than 
AA collisions in order to test the predictions of the color glass condensate, 
because the details of the system at early times are carried out to the final state 
almost unaltered. One should therefore be able to perform more direct tests of 
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the CGC ideas in the context of dA collisions. 



4.4.2 Anisotropy 

Another interesting quantity to look at is the so-called "elliptic flow" , which 
signals the existence of a significant pressure in the transverse direction [56], 
that converts the original spatial anisotropy into an anisotropy in momentum 
space. 

This has been estimated by Krasnitz, Nara and Venugopalan [57] from the 
numerical solution of the classical Yang-Mills equations. The growth of the 
anisotropy with time is similar to what happens in hydrodynamical models, [57] 
and the full v 2 is attained after a time of the order of the size R of the system. 
The V2 obtained at late times in the color glass condensate framework is plotted 
in figure 7. One can see that although it is not in quantitative agreement with 



Figure 7: The centrality dependence of V2 from the classical gluon production 
(green band). The black filled dots are STAR data. 

the data, it has the correct order of magnitude and reproduces qualitatively the 
dependence on the multiplicity (i.e. on the centrality). However, the dependence 
of V2 on the transverse momentum is in clear disagreement with the RHIC data: 
the color glass condensate predicts that V2 has a maximum at small momentum 
and then falls to zero at large momentum, while the measured V2 increases 
monotonously up p± ~ 2 GeV and remains almost constant at higher momenta. 
Note that the description of the system in terms of classical fields breaks down 
for times that are large compared to Qj 1 , that is for times needed to establish 
the V2- It is nevertheless instructive to see that this model, which does not rely 
on the assumption of thermalization, produces a V2 of comparable magnitude 
to the «2 obtained in hydrodynamical expansion over similar time-scales. 
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4.5 Ratios of spectra: Raa, RdA 



The comparison of the p^ spectra in AA or dA to the corresponding ones in pp 
collisions is usually done through the following ratio: 



<//v 
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dyd' J p 1 
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dyd 2 i 
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where N co \\ is the number of binary collisions. For the production of high p± 
particles, expected to scale like the number of binary collisions, these ratios 
should be unity. These ratios have been measured at RHIC both for AA and 
dA collisions. 
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Figure 8: The ratio RdA measured at RHIC by the BRAHMS collaboration at 
several rapidities (positive rapidities correspond to the fragmentation region of 
the deuteron) [58]. 

In A A collisions at central rapidity (rj = 0), the measured Raa is much 
smaller than one (around 0.2) in central collisions, even at pj^'s as high as 
10 GeV, and approaches one from below as one goes to more and more peripheral 
collisions. However, at the values of x probed by gluon production at mid- 
rapidity (x ~ 10~ 2 for p^ ~ 1 GeV), where one expects the distribution M^ [p] 
to be given by the McLerran-Venugopalan model (eq. (8)), the Raa predicted 
in the color glass condensate at a time of the order of Qj 1 has a maximum 
larger than 1 at small p± and then goes to one from above as p± increases 
[59]. This maximum of Raa in the MV model is interpreted as a manifestation 
of the multiple (elastic) rescatterings of the produced gluon, which tend to 
redistribute the p± spectrum by depleting the small pj_'s and enhancing the 
high p±'s (Cronin effect). The observed suppression of Raa therefore requires 
final state effects, and is naturally interpreted as energy loss ([60-64], see also 
[65] for a recent review) of the high p± partons as they go through the dense 
medium formed in A A collisions. 

At much smaller values of x, which can be probed by looking at particle 
production at large rapidity, one expects the distribution W Xo [p] to reach the 
asymptotic form of eq. (7). It has been found that with such a distribution 
of color sources, the color glass condensate leads to a suppression of the ratios 
Raa or RdA [66] (see the plots of figure 9, taken from [67]). However, in order 
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to disentangle this suppression which is an initial state effect (because it comes 
from the "wave function" of the projectiles) from the energy loss, one needs an 
experiment in which one has no significant final state effects, like dA collisions. 
The ratio RdA measured by the BRAHMS experiment is shown in figure 8 [58]. 
One can see at r\ = a very different behavior than for A A collisions: the ratio 
RdA has a maximum above one and seems to go to one from above at high 
p±. As one increases the rapidity of the observed particles, one can see the 
ratio RdA drop very fast and eventually become consistently smaller than one. 
Such a behavior with rapidity was predicted in [67] (see the figure 9), where 
the evolution of the ratio RdA was evaluated by evolving the unintcgratcd gluon 
distribution in a nucleus using the Balitsky-Kovchegov equation, starting from 
the McLerran-Venugopalan model at large x. Similar qualitative results were 




Figure 9: Evolution of the ratios Raa and RdA with rapidity in the color 
glass condensate framework. The initial condition is the McLerran-Venugopalan 
model. Note that the calculation of Raa assumed k± -factorization. 

also obtained in [68] , in [69] within a toy model for the gluon distribution inspired 
from the color glass condensate, and in [40] where the two distributions of frozen 
sources given in eqs. (7) and (8) were used in order to compute the ratio RdA- 
The onset of the suppression of the ratio RdA as one increases the rapidity was 
studied analytically in [70] , where it was shown that the rapid suppression seems 
to be due to the fact that the gluon distribution evolves faster in a proton than 
it does in a nucleus, because the nucleus reaches the saturation regime earlier. 
Recently, a more quantitative study of this effect has been performed in [71]. So 
far, the color glass condensate is the only framework in which one reproduces, 
at least qualitatively, both the the Cronin effect at central rapidities and its 
suppression at forward rapidities. 
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5 Discussion and outlook 



Considerable theoretical progress has been made over the last ten years in un- 
derstanding the physics that govern the wave function of a hadron at high 
energy. Not only has one acquired a fairly intuitive picture that is universally 
applicable to all hadrons at high energy, but an operational framework has been 
developed in which many phenomena can be described quantitatively One is 
now in principle able to study and make predictions in the entire x, k± plane 
(with the exception of the truly non-perturbative region k± < A QCD ), with non- 
perturbative physics entering only as an initial condition for otherwise known 
evolution equations. 

With the realization that the characteristic scale that governs high energy 
hadronic interactions, the saturation momentum Q s , is enhanced by the size 
of the projectiles in collisions involving nuclei, RHIC can be viewed as a place 
of choice in order to test these ideas and confront them to data. And indeed 
the phenomenology based on saturation physics has been quite successful in 
describing what is observed at RHIC. Whenever relevant comparisons with data 
can be made (global observables in AA collisions, spectra in dA collisions), the 
tests are successful. In cases of marked disagreement, like with the elliptic flow, 
one understands that the discrepancy comes from using the color glass picture 
beyond its domain of applicability. 

One may ask whether the predictions which have been tested so far are 
distinctive features of the color glass condensate that cannot be reasonably 
explained by other models. Naturally, the more precise the question, the more 
discriminatory the answer. Because in AA collisions the final state interactions 
are important, only global observables are preserved from early times to the 
final state. Therefore AA collisions are in general not ideal for direct tests of 
the color glass condensate (many global features of AA collisions are present 
in any reasonable model of hadronic interactions, and thus not characteristic 
of the color glass condensate). Collisions involving a small projectile, like dA 
collisions, where effects of final state interactions can be minimized, at least in 
some kinematical domains, are much better suited. And indeed the suppression 
of the ratio RdA at forward rapidities could turn out to be such a distinctive 
measurement, since no other model has been able so far to provide a natural 
explanation of what is observed. This is also one instance where non trivial 
quantum evolution could be playing a significant role. What is missing for 
this to become a real test is a detailed calculation of this effect that would go 
beyond the many qualitative approaches performed so far. Clearly, doing such 
calculations, as well as getting more precise data, is of utmost importance. 

From comparisons with data, one has also learned that final state inter- 
actions are generally very important in AA collisions, and many observables 
do not simply reflect the properties of the partonic system produced in the 
early stages. The color glass condensate only provides the initial condition for 
a subsequent evolution of the system leading possibly to the formation of a 
quark-gluon plasma. Understanding whether and how thermalization happens 
presents interesting theoretical challenges (see [72-76] for some recent works on 
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the subject) and is of utmost importance for giving a solid theoretical basis to 
present descriptions of AA collisions. 

As a final remark, oriented towards the future, one may note that many of 
the phenomena uncovered at RHIC, should become more clearly visible at the 
LHC. There, with center of mass energies of 5.5 TeV for AA collisions, the typical 
value of x at mid-rapidity will be about 5.1CT 4 and values as small as 10~ 5 
could be reached at forward rapidities. This corresponds roughly to values of Q 2 S 
ranging from 5 GeV 2 to 14 GcV 2 . Such large values of the saturation momentum 
make the coupling constant a s smaller than at RHIC, giving firmer grounds to 
the weak coupling expansion used in the color glass condensate framework. A 
corollary of this is that the gluon occupation number in the saturation region 
will be larger, making the classical description also more justified. 
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